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I. In a recent Letter Yang and Bodek jjj presented 
results of a new analysis of proton and deuteron structure 
functions in which the free neutron structure function, 
F£, was extracted at large x. Knowledge of Ff is crucial 
for determining the neutron/proton structure function 
ratio, whose x — ► 1 limit is sensitive to mechanisms of 
SU(6) spin-flavor symmetry breaking, and provides one 
of the fundamental tests of the x dependence of parton 
distributions in perturbative QCD. 

Relating nuclear structure functions to those of free nu- 
cleons is, however, not straightforward because at large 
x nuclear effects become quite sizable. In particular, 
omitting nuclear binding or off-shell corrections can in- 
troduce errors of up to 50% [^],[| in F^/F^ already at 
x ~ 0.75. Rather than follow the conventional proce- 
dure of subtracting Fermi motion and binding effects in 
the deuteron via standard two-body wave functions [Q, 
Yang and Bodek instead extract FT 1 using "a model pro- 
posed by Frankfurt and Strikman H , in which all binding 
effects in the deuteron and heavy nuclear targets are as- 
sumed to scale with the nuclear density" Q. Here we 
point out why this approach is ill-defined for light nuclei, 
and introduces a large theoretical bias into the extraction 
of F£ at large x. 

For heavy nuclei the nuclear EMC effect is observed to 
scale with the nuclear density, pa (|-@] 
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where R A = PAjF^ and p A = SA/^ttR^), with R\ = 
(5/3)(r 2 ) and (r 2 ) 1 ' 2 is the nuclear r.m.s. radius. As 
suming that an analog of Eq.(Q) holds also for F^/F^ 
(F-f = Ff + F%) one finds 
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This expression was derived by Frankfurt and Strikman 
in Ref. [pi, where the denominator was further approx- 
imated by pa — Pd Ra ~ PA — Pd- It was used in the 
analysis of the SLAC data jfj), and also referred to by 
Yang and Bodek , although the explicit formulas used 
in that analysis are not given. We have checked the nu- 
merical values for F^/F^f in [jjj, and they agree with 
the result one would obtain from Eq.(||) if the densities 
quoted in Q are calculated in terms of charge radii ||. 

Frankfurt and Strikman point out || that from the 
above expression for F% /F% one can extract the free neu- 
tron structure function from empirical EMC ratios and 



the nuclear densities. With the numerical values for pa 
quoted in |J, one finds then that the EMC effect in d is 
about 25% as large as in 56 Fe |j]||, and has the same x 
dependence. 

While the correlation of EMC ratios with nuclear den- 
sities is empirical for heavy nuclei, application of Eq.(jl]) 
to light nuclei, A < 4, for which EMC effect has not yet 
been determined, is fraught with ambiguities in defining 
physically meaningful nuclear densities for few-body nu- 
clei. Firstly, the relevant density in Eq. (|lj) is the nuclear 
matter density, while in practice pa is calculated from 
the charge radius [|l]J(| — for heavy nuclei the difference 
is negligible, but for light nuclei it can be significant. 
Secondly, treating the deuteron as a system with radius 
(r 2 ) 1 / 2 ss 2 fm means that one includes both nucleons 
in the average density felt by one of them, even though 
one nucleon obviously cannot influence its own structure. 
Therefore what one should consider is the probability of 
one nucleon overlapping with the other, which is simply 
the deuteron wave function at the origin. This has zero 
weight, however, so the only sensible definition of mean 
density for the deuteron is zero. Strictly speaking, the 
nuclear density extrapolation then predicts no nuclear 
EMC effect in the deuteron. 

In Ref. ^ Frankfurt and Strikman argue that for heavy 
nuclei the average potential energy is proportional to the 
average nuclear density, and hence for x below 0.5-0.6 
(where nuclear Fermi motion is not overwhelming) the 
nuclear EMC effect should scale with average nuclear 
density. If one applies the idea from heavy nuclei (where 
the assumption is known empirically to be reasonable) 
to the deuteron, one finds that the EMC effect in d is 
(Fi/Fi 1 - 1) = 0.25 (Ff e /F£ - I). For light nuclei 
(A = 2, 3) no justification for this assumption is provided, 
however, and for x > 0.6, where nuclear Fermi motion ef- 
fects become large, Frankfurt and Strikman caution that 
this estimate is only a qualitative one 0. 

In a reply to our original Comment, Yang and Bodek 
state that "although the notion of nuclear density for the 
deuteron may not be very well defined, the value of the 
nuclear density for deuterium that was used in the SLAC 
fit yields a similar correction for nuclear binding in the 
deuteron as the estimate by Frankfurt and Strikman" 
Q. As explained above, not only is the notion of nuclear 
density for the deuteron not very well defined, it is not 
defined at all. 

Moreover, agreement with other calculations [|6|]l0|] for 
the magnitude of the deuteron EMC effect does not 
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provide a posteriori justification for using an ill-defined 
quantity like average nuclear density for the deuteron. 
One would never think of using a density extrapolation to 
extract the neutron's electromagnetic form factors from 
quasi-elastic scattering on the deuteron or 3 He, for ex- 
ample, and there is no reason to believe this method is 
any more reasonable for structure functions. 

II. The size of the EMC effect in the deuteron cannot 
be tested directly in any inclusive deep-inelastic scatter- 
ing experiment on the deuteron, as it requires knowledge 
of F£, which itself must be extracted from deuteron data. 
If, on the other hand, the EMC effect scales with nu- 
clear density even for the deuteron, as assumed in [jl],^), 
it must also scale with pa for all A > 2. In particu- 
lar, it must predict the size of the EMC effect in 3-body 
nuclei. In fact, for A = 3 the nuclear density extrap- 
olation makes quite a dramatic prediction: since the 3- 
body nuclear densities calculated from the charge radii 
are p 3He = 0.049 fm" 3 and ps H = 0.068 fm~ 3 , the EMC 
effect in 3 H is 40% larger than that in 3 He. This is 
to be compared with standard many-body calculations 
in terms of Faddeev wave functions [|ll| which predict 
< 10% difference between the EMC effects in A = 3 mir- 
ror nuclei, see Fig.l. The A = 3 system presents therefore 
an ideal case for testing the scaling of the nuclear EMC 
effect to small A. 



Density extrapolation 




FIG. 1. Ratio of EMC effects in 3 He and 3 H using standard 
many body (Faddeev) wave functions, and the nuclear density 
extrapolation. 



III. In Ref. H] , Yang and Bodek compare their F% / 
ratio with the relativistic quark model calculation of Ref. 
fiof (which we denote by "MST" , and which is referred 
to in jlj] as the "Melnitchouk-Thomas theory"), while 
attributing to this the empirical extraction of F^ which 
was actually carried out in Ref. M (which we denote by 
"MT"). They fail to appreciate the distinction between 
the MST investigation of relativistic and off-shell effects 
within a model and the MT analysis of data, so a brief 
explanation is needed to clarify the confusion. 

In the MST model, the objective was to investigate 
the extent to which relativity and nucleon off-mass-shell 



effects (p 2 ^ M 2 ) violate the familiar convolution ap- 
proximation for nuclear structure functions. To test the 
significance of these effects, a microscopic model of nu- 
cleon structure functions had to be constructed (in the 
absence of data on off-shell nucleon structure functions!), 
in order to consistently combine this with a relativistic 
deuteron wave function [l2| ] in a covariant calculation. 
Both the proton and neutron structure functions were 
parameterized in terms of a simple model for the nucleon- 
quark-diquark vertex functions, and, in the context of 
this model, the theoretical F^/F^ ratio was constructed. 

No attempt to extract F% from the calculated ratio was 
made in the MST model. Indeed, if any input neutron 
structure function is used to construct the theoretical 
F2/F2 ratio, clearly the calculated ratio cannot then be 
applied to draw conclusions about the extracted F% - 
different neutron input would result in a different EMC 
ratio, and the argument would be cyclic. Therefore com- 
paring the ratio calculated in the MST model with a ratio 
used in an empirical extraction of F£ such as in Ref. Q| 
is rather meaningless. 

The least ambiguous procedure is to deconvolute the 
smeared neutron structure function by iterating the in- 
version procedure until a convergent, self-consistent so- 
lution is obtained. This was done in the MT analysis, in 
which the only theoretical input was the deuteron wave 
function. 

In Ref. §, Yang and Bodek criticize the calculated 
F^/F^ ratio in the MST modelQ on the grounds that 
its small- a; behavior is opposite to that in the EMC ra- 
tio for iron. They also assert that the small- a; behavior 
reflects some sort of violation of the energy-momentum 
sum rule in the MST model. The history of this discus- 
sion is rather long and illustrious, and for the background 
we refer the reader to any standard review of the EMC ef- 
fect (see e.g. Refs. ]7|,|l3[ and references therein). Both of 
these points are actually diversions and irrelevant to our 
main comment, however, since they were raised as ap- 
parent justifications for using the density extrapolation 
model, we feel a need to address them. 

Firstly, the behavior of structure functions at x ~ 0.2 
is clearly outside the region of large x where the nuclear 
effects which we are discussing are relevant. The focus of 
the MST investigation was specifically nuclear effects in 
the valence quark region at large x. A realistic descrip- 
tion of the sea requires much more sophistication than 
can be reasonably demanded of any simple valence quark 
model. Secondly, in any realistic model of nuclear struc- 
ture which properly accounts for nuclear binding, nucle- 



1 Note that this is incorrectly attributed here to high- 
lighting the apparent confusion in between the MST 
model and the MT analysis. The curve labeled "Melnitchouk- 
Thomas model" in Fig.l of jl]] is also incorrectly referenced 
to H rather than to the MST model @. 
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ons alone obviously cannot carry all of the momentum 
of the nucleus. An analog of demanding that nucleons 
saturate the momentum sum rule would be to demand 
that valence quarks alone carry all of the nucleon's mo- 
mentum and gluons carry none — in contradiction with 
experiment. This is the content of the assertion || that 
the energy-momentum sum rule is not satisfied in the 
MST model. 

The point is that we have never advocated using the 
MST model to extract the neutron structure function, 
contrary to the claim made in As explained above, 
the least ambiguous F£ extraction requires employing 
the iterative procedure for unsmearing the effects of 
the deuteron wave function, as outlined in ||||, using 
deuteron and proton data as input. However, even if 
these points were valid, they would still not be relevant 
to the issue of F^ 1 extraction as they refer specifically to 
the MST model and not to the F 2 ™ extraction in the MT 
data analysis. 

Yang and Bodek in addition claim P] that the MST 
model has not been applied to A > 4 nuclei. Since the 
MST model was relativistic (namely, it included terms 
beyond order v/c in the deuteron wave function), ap- 
plying it to other nuclei is straightforward once rela- 
tivistic wave functions are known. To date only non- 
relativistic wave functions have been calculated for nu- 
clei with A > 4. However, since the MST model has a 
well-defined non-relativistic limit |Q (namely, omitting 
deuteron P-waves, and dropping terms of order v 2 /c 2 ), 
it smoothly matches onto previous non-relativistic calcu- 
lations of the nuclear EMC effect for A > 4 So 
again, contrary to the assertion by Yang and Bodek |j| , 
the MST model has indeed been tested for all nuclei and 
for all cases where wave functions exist [fL5[ . 

IV. To summarize, we have demonstrated that the nu- 
clear density model of the EMC effect, extrapolated well 
beyond its region of validity to the case of the deuteron, 
is a completely unreliable method of extracting the neu- 
tron structure function at large x, and introduces a large 
theoretical bias into the extraction procedure. Although 
a nuclear density fit may be quite useful for heavy nuclei 
where data exist ||, its extrapolation to A < 3 where 
there are no data is highly speculative, as illustrated by 
the difficulty in defining physically meaningful densities 
for few-body systems. 

We believe this issue can only be resolved by measur- 
ing the nuclear EMC effect in the deuteron and in A = 3 
nuclei. One proposal |l6| would be to simultaneously 
measure the structure functions of 3 He and 3 H, extract- 
ing the Fg /F% ratio through the cancellation of nuclear 
effects, which would indirectly determine the EMC effect 
in the deuteron. Other solutions would be to reconstruct 
F£ from the d/u ratio extracted from parity-violating ep 
scattering p7|, completely free of nuclear effects, semi- 
inclusive 7T^~ production from an X H target JTq] , or from 
various charged-current reactions [fL9| . The above discus- 



sion presents a strong case for performing these experi- 
ments as soon as possible. 

We would like to acknowledge stimulating discussions 
with G. Petratos concerning the nuclear corrections in 
^4 = 3 nuclei. This work was supported by the Australian 
Research Council. 
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